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We explore the possibility to construct higher-twist parton distributions in a nucleon at some 
low reference scale from convolution integrals of the light-cone wave functions (WFs). To this end 
we introduce simple models for the four-particle nucleon WFs involving three valence quarks and 
a gluon with total orbital momentum zero, and estimate their normalization (WF at the origin) 
using QCD sum rules. We demonstrate that these WFs provide one with a reasonable description 
of both polarized and unpolarized gluon parton densities at large values of Bjorken variable x > 0.5. 
Twist-three parton distributions are then constructed as convolution integrals of qqqg and usual 
three-quark WFs. The cases of the polarized structure function g2{x,Q 2 ) and single transverse 
spin asymmetries are considered in detail. We find that the so-called gluon-pole contribution to 
twist-three distributions relevant for single spin asymmetry vanishes in this model, but is generated 
perturbatively at higher scales by the evolution, in the spirit of GRV parton distributions. 
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I. INTRODUCTION 

Higher-twist parton distributions are conceptually very 
interesting as they go beyond the simple parton model de- 
scription and allow one to quantify correlations between 
the partons. Unfortunately, they prove to be very elu- 
sive. Despite considerable efforts, very little is known 
even about the simplest, twist-three distributions which 
contribute, e.g. structure function gi(x, Q 2 ) in the polar- 
ized deep-inelastic scattering [THS] and transverse single 
spin asymmetries (SSAs) in the collinear factorization 
approach [TOrEO] . 

One general reason for this is that the structure of 
higher twist parton distributions is much more compli- 
cated compared to the leading twist: they are functions 
of two and more parton momentum fractions. The usual 
strategy to extract parton distributions from experimen- 
tal data has been to assume a certain functional form 
with a few adjustable parameters at a reference scale, and 
find the parameters by making global fits to the available 
data. This is a standard approach which works quite well 
for the leading twist. Unfortunately, it does not work for 
higher twist (or, at least, has not been applied systemat- 
ically) because there is no physical intuition on how such 
distributions may look like. Also the asymptotic behav- 
ior of higher-twist distributions both at small and large 
x is poorly understood. Hence it is very hard to guess an 
adequate parametrization. 

In this work we make a step in this direction. Recall 
that the case of higher-twist parton distributions is not 
unique in that they are functions of several kinematic 
variables: in studies of generalized parton distributions 
(GPDs) or "unintegrated" transverse-momentum depen- 
dent distributions (TMDs) the same complication arises. 
In both cases, representations in terms of overlap inte- 
grals of light-cone wave functions have been extremely 



useful for developing the underlying physics picture and 
provide one with a good basis for theoretical modelling. 
In what follows we try to follow the same path for the 
construction of higher-twist distributions as overlap in- 
tegrals between Fock states with the minimum (valence) 
and next-to-minimum (one extra gluon) parton content. 

In order to keep the model as simple as possible, in this 
work we restrict ourselves to contributions of the states 
with total zero angular momentum. We overtake the 
expressions for three-quark wave functions from Ref. [21] 
which have been shown [2H [22] to provide one with a 
good description for quark parton densities at large x and 
the nucleon magnetic form factor. The new contribution 
of this paper is to include into consideration the Fock 
states with one additional gluon which were considered 
in [21] on a qualitative level. We find that there exist 
three independent qqqg wave functions with zero orbital 
momentum. Our analysis of their symmetry properties 
does not agree with earlier results [23]- We calculate the 
normalization of these new wave functions using QCD 
sum rule approach and construct explicit models by the 
requirement that their light-cone limit (zero transverse 
separation) reproduces the nucleon twist-4 distribution 
amplitudes introduced in Ref. [24] . 

Having specified the wave functions, we calculate 
the quark and gluon polarized and unpolarized par- 
ton distributions and find agreement with the existing 
parametrizations at large x without any fine-tuning of the 
parameters. Encouraged by this, we construct the twist- 
three correlation function involving a quark, antiquark 
and gluon fields which is relevant for the structure func- 
tion g2{x, Q 2 ) and single spin asymmetries. In our model 
this correlation function vanishes at the boundaries of 
parton regions where one of the momentum fractions goes 
to zero, but non-zero values are obtained at higher scales 
perturbatively through the QCD evolution. This phe- 
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nomenon is in full analogy to the generation of a large 
gluon parton distribution at small x starting from the 
"valence" -like ansatz in the GRV approach [25,(25]. Such, 
radiatively generated, soft-gluon pole and soft-fermion 
pole contributions to the spin asymmetries are calculated 
and compared to the existing parametrizations. The sign 
of radiatively generated soft pole terms as well as the sign 
of the twist-three contribution to the structure function 
92{x,Q 2 ) at large x are largely model- independent pre- 
dictions of our approach; these signs turn out to be in 
agreement with the data in all cases. Finally, we discuss 
possible generalizations of our simple model that may 
provide one with usable parametrizations for the phe- 
nomenological analysis. 



II. LIGHT-CONE COORDINATES 

For an arbitrary four- vector we define the light-cone 
coordinates as 

a+ = -L( a ° + a 3 ), a_ = -L( a ° - a 3 ) , 



a = a 1 + ia 2 



a = a 1 — ia 2 



(1) 



so that the matrix a = a^a^, where cr^ = (I, (?) takes the 
form 



'2a,- -a 
—a \p2a^ 



(2) 



In what follows we use Weyl representation for the 
7— matrices 







a' 



I ¥ I ' 7 V -cr* 



5 _ „,0„,1„,2„,3 / — I 



7 = 7777 







(3) 



and the two-component notation for Dirac spinors 



i'a) = (Ql 



(4) 



The two independent light-like vectors 



1 



1 



= -=(1,0,0,-1), ft" = -=(1,0,0,1), (5) 



V2 



V2' 



n 2 = n 2 = 0, nn = 1 can be parametrized in terms of the 
two auxiliary Weyl spinors: 



^ao X a Xa ) ^ad — f^af^o 



(6) 



where 



K = 2 1 / 4 , ^ = 2 1 / 4 (J) , 

a, = , fe =2v*(;y (7 ) 



We accept the following rules for raising and lowering the 
spinor indices (cf. Ref. [24] ) 



A™ 



c a0 



Xn 



X 



£/3a 



x c 



A /3 e 1 



Xr 



(8) 

where the antisymmetric Levi-Civita tensor is defined as 

e 12 = e 12 = -e^ = -e i2 = l. 

The auxiliary spinors A and /x are normalized as 

(vX) = n a X a = -(Xfi) = -V2, 

(/2A) =^A d = -(A/l) = +V2 (9) 

and serve to specify "plus" and "minus" components of 
the fields. We define 



x+ 



x&x a : x- 



Xafi a , 



(10) 



so that each two-component spinor can be decomposed 



as 



QiX)lpa = Aq 1p- 

(A/2)x d = A d x- 



P- a x+- 



(11) 



In the same notation the light-cone decomposition of a 
vector (e.g. gluon) field takes the form 

A A 

A a a = ^4— X a Xa + A+ l-l a p,a H 7= X a jla H T? ^aXa ■ 

v 2 v 2 

(12) 

The " plus" spinor fields ip+ , x+ an d transverse gluon 
fields A, A are assumed to be the dynamical fields in the 
light-cone quantization framework. The "minus" fields 
V>_ , X— i A— can be expressed in terms of the dynamical 
ones with the help of equations of motion (EOM) whereas 
A + = due to the gauge fixing condition. 

The plus quark fields have the following canonical ex- 
pansion 

dp + d 2 p± 



^ (2tt)3 
dp + d 2 p± 
v/2pT(2^) 



e- tpx h(p) 
e- ipx b t (p) 



5 +lpX d\(p) 

z +ipx d\(p) 



(13) 

where &f(n, ^t(^) are l ne annihilation operators of quark 
and antiquark of positive (negative) helicity, respectively. 
They obey the standard anticommutation relations 



{b x (p),bUp')} = {dx( P ) > dUp')} = 
= 2p + (2ir) 3 6 Xi y6( P+ -p' + )5 2 (p A 



P±) 



(14) 

The similar expansion for the dynamical transversely po- 
larized gluon fields A and A reads 



t/ x f dk+ d 2 k± . 

^ = ^/4(#^ 

, , x f dk + d 2 k± , 

^ = ^74(2^^ 



*a t (fc) 



,,+ifcx t 



(k) 



—ik 



x ai {k) + e +lkx a\{k) 



(15) 
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Here and below A — ^2, a t a A a etc. where t a are the 
usual SU(3) generators in the fundamental representa- 
tion, normalized as tr(t a t b ) = ^ ab - The creation and 
annihilation operators obey the commutation relation 

~a b x(p),(4(p'))i = 

= 2p + (2Tr) 3 6 x , y 5 bb 'S(p + -p' + )S 2 (p ± -p' ± ). (16) 

Finally, the gluon strength tensor F^ v and its dual F^ v 
can be decomposed as 



F 



a/3,a/3 



IF 



(17) 



Here / Q| a and f^i are chiral and antichiral symmetric 
tensors, f a p = fp a , f = /*, which belong to (1,0) and 
(0, 1) representations of the Lorenz group, respectively. 
Their "good components" are defined as 



/++ — A Q A /3 / Q , ( g , 
In the light-cone gauge 

f ++ = -d+A, 



/++ = A^/~ 



a/3 



f ++ = -d + A, 



(18) 



(19) 



where 9+ = n^d^ = d/dx-, so that they can readily be 
expanded in contributions of annihilation and creation 
operators using Eq. (15). 



As mentioned above, "minus" field components can be 
expressed in terms of the dynamical fields using QCD 
equations of motion. 



III. NUCLEON LIGHT-CONE WAVE 
FUNCTIONS 

A. Definitions and symmetry properties 

The light-cone wave functions (LCWFs) are defined 
as probability amplitudes of the corresponding parton 
states which build up the proton with a given helicity. 
They depend on parton longitudinal momentum fractions 
Xi, transverse momenta fcj^, and on parton helicities. 
LCWFs are usually thought of as solutions of the eigen- 
value problem for the light-cone quantized QCD Hamilto- 
nian [37J HE] > although this construction is far from being 
complete. 

Throughout this work we adopt some definitions and 
partially also the notation from Ref. [21] . In particular 
we use a shorthand notation for the 7V-parton differential 
phase space 

A' 



[dx] N = Y[dXiS(l- j^Xj) , 



N 



[dk±]N = (16*V-i II * k ±* 52 (E (20) 



and 



[VX] 



N 



[dx}N[dk±] 



±\N ■ 



(21) 



The valence three-quark state with zero angular mo- 
mentum is the simplest one. It can be described in terms 
of the single LCWF [2TJ 130] 



' uud 



e ijk 



(22) 



(«| r (l)«^(2)4 t (3) - 4 f (l)4(2)4 f (3)) |0) . 



Here and below the argument of the field w]^.(l) etc. 
refers to the collection of its arguments that arc not 
shown explicitly, i.e. uL(l) = vlAxi, k± t i). The (real) 



function ^^ 3 (X) depends on momentum fractions Xi 
and transverse momenta of all partons. 

Models for *W (X) of various degree of sophistication 
have been was considered in different context in a large 
number of papers see e.g. Refs. [2lTf2"3l I29H31"] . In this 
work we adopt the simplest ansatz [21] 



(o) _ 

123 



1 

'WE 



(x 1 ,x 2 ,x 3 )Q 3 (a 3 ,x i ,k 1 _ i ) . (23) 



The transverse momentum dependence is contained in 
the function 0,^ 



il N {a N ,x 2 ,kj_i) = —— — exp 



X\X 2 



.x N 



'■N 



(24) 



which is normalized such that 



[d 2 k±] N n N (a N ,x i ,kj_ l ) = 1, 
[d 2 k±] N n 2 N (a N , Xi, k±i) 



Pn 



Xi 



■XN 



(25) 



where 



Pn = (87r a 



N) 



and <j>(Xi) is related to the leading-twist-3 nucleon distri- 
bution amplitude (see the next Section). The parame- 
ter 03 determines the spread of the wave function in the 
transverse plane and e.g. the average quark transverse 
momentum. 

The general classification of Fock states involving an 
additional gluon was given in Ref. |23j . Unfortunately, 
we do not agree with the analysis in [53] of the symmetry 
properties of the corresponding LCWFs. 

As in the three quark case, we restrict ourselves to the 
states with zero total orbital angular momentum, L z = 0. 
There are two possibilities [23] : either the quark helicities 
sum up to X uu d = 3/2 and the gluon has opposite helicity 
to that of the proton, X g = — 1, or, alternatively, X uu d — 
— 1/2 and \ g = +1. We begin with the first case. 
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The starting observation is that the SU(3) generators 
obey the following identity 



definite parity under cyclic permutations of the first three 
(quark) arguments 123 — > 231: 



e ijl tf k + e ilk tfj + e ljk t% = 



(26) 



As a consequence, there exists only one possibility to 
form a colorless state (up to equivalent redefinitions) 



= e ijk J[DX] 4 ^ 1234 (X) 
x 5 f(4)[t°u t (l)]ln]. t (2)4 t (3)|0). (27) 

Note that [i a w^(l)]| = uL.(l)i^. Symmetry properties 

of the LCWF *f 234 are determined by the requirement 
that the nucleon has isospin 1/2. Since I3 = 1/2 is fixed 
by the quark flavor content, the 7 = 1/2 requirement 
is equivalent to the simpler condition that the state is 
annihilated by the isospin step-up operator 



I+\P, 



0. 



The action of I + amo unts to the replacement of quark 



flavors d — > u in (27), I 



S/6S. 



Projecting the 
resulting state onto (0\gf {A')uk'^(i')uji^(2')uki^(V) and 
collecting the terms in the two independent color struc- 
tures (cf. Eq. (26)) one finds two constraints: 



*i234 + *1324 



2134 



2314 



^3124 
llf J- 

^3124 



*3214 



3214 



= 0, 

= 0. 



(28) 



Since the second equation can be obtained from the first 
one by renaming 1 O 2, only one of them is indepen- 
dent. In order to solve this constraint it is convenient to 
represent the function as sum of contributions with 



^1234 



^1234 



^1234 



*i324 > 



(29) 



such that 



*1234 



= W 



^,0 
2314 ! 



T 1234 — e ^2314' 



One easily finds that an arbitrary function ^1234 1S 



*1234 



"^1324- Thus the most general solution to the 



solution of Eq. (28) whereas one has to require that 



isospin constraint can be written as 



*1234 



v 1234 



*1234 



*1324 i 



(30) 



where and + are arbitrary functions with the 

specified symmetry under cyclic permutations. 

Our result does not agree with the conclusion of [23 



that the function ^\ 2 34 (^uudg m notations of Ref. [23]) is 
antisymmetric with respect to permutation of the second 
and third arguments, which is a much stronger condition. 
In fact any function which is antisymmetric in 2 f> 3 
can indeed be written in the form (30). However, e.g. a 



totally symmetric function in the quark arguments is also 
allowed. The reason why this does not contradict isospin 
counting is that the corresponding state is annihilated by 
I + thanks to the color identity ( 26 ) . We note in passing 



that the SU (3)-color generators in the definitions given 
in [53] must be transposed, tf t , — > t° v 

The second case, a gluon with positive helicity, can be 
treated similarly. There exist two independent LCWFs 
which can be defined as 



\Pi +)uudgt 



ijk 



Mif^W [* a «4.(l)]j(«J t (2)4 4 .(3) - 4 t (2)«^(3)) 5 -' t (4) 

(JT)«t ; (i)([t- u+ (2)]t4 1 .(3) - [i-d ; (2)]t«t r (3))^' t (4)}|0) . 



T(2) , 

1234*. 



(31) 



1~ ( 1 1 

The functions ^x234 anc ^ ^1234 have no symmetry con- 
straints. This result also does not agree with [23 . 

In what follows we accept the following ansatz for the 
quark-gluon LCWFs 



t(2) 



*1234 



*t(l) _. 
*1234 — 

*t(2) 
*1234 — 



1 



/2S4 
1 

1 

/2X4 



(f> g (xi , x 2 , x 3 , ar 4 ) fi 4 (a^ , Xi , fcj_i ) , 



if> g 2) (xi , x 2 , x 3 , X4) fi 4 (a* ,Xi,k±i). (32) 



The function Q4 is defined in Eq. ( 24 ) and the momentum 



(1 2) 

fraction distributions <f) g (xi), ip g ' (a;,) are related to the 
next-to-leading twist-4 nucleon distribution amplitudes 
as discussed in the next section. For simplicity we choose 
the same parameter al determining the spread of both 
wave functions v]/t(i) an d \j/t( 2 ) j n the transverse plane. 
This restriction can be relaxed. 



B. Relation to nucleon distribution amplitudes 

Nucleon distribution amplitudes (DAs) are defined as 
LCWFs with all constituents at small transverse separa- 
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tions, schematically [3"U] 

n\k±\<H 

(t)(xi,n)~l [dk±] N ® N (x i ,k_L ! i) . (33) 

As always in a field theory, taking an asymptotic limit 
(here vanishing transverse distance) produces diver- 
gences that have to be regularized. Hence DAs are scale- 
dependent objects which only include contributions of 
small transverse momenta, less that the cutoff. 

The exponential ansatz for the the transverse momen- 



tum dependence of the LCWFs ( 23 1 , (32 1 implicitly as- 



sumes that contributions of hard gluon exchanges *~ l/kj_ 
are subtracted as well, so that it is natural to identify in- 
tegrals of the LCWFs over transverse momenta with the 
corresponding DAs at a certain low normalization scale. 
The advantage of of imposing this condition is that nu- 
cleon DAs allow for a different and more rigorous defi- 
nition in terms of matrix elements of nonlocal light-ray 
operators. Their moments can be studied using Wilson 
operator product expansion (OPE) and estimated using 
QCD sum rules and/or lattice calculations. The identifi- 
cation of the integrals of the LCWFs with (dimcnsionally 
regularized) DAs can be viewed as the choice of a specific 
renormalization (factorization) scheme. 

To begin with, consider the leading twist-three nucleon 
DA which is defined by the matrix element [35] 

(0|e y ' fc (ul T {z 1 n)C^(z 2 n)^dl(z 3 n)\p) = 
= -\ P+ iN\p) /[dz] 3 e-*+5>*** 3 (:B<)> (34) 



where N(p) is the nucleon Dirac spinor, p 2 — m 2 N , 
N^(p) = |(1 + 75) AT (p) and C is the charge conjuga- 
tion matrix. Going over to the two-dimensional spinor 
notation and using the explicit expression for the C- 
matrix in Weyl representation |33) 



C = z 7 2 7° = ( 









e a0 



(35) 



this definition can be rewritten equivalently as 

(0\e i ^(z 1 )u^(z 2 )d^(z 3 )\p, +> = 



V2 



.3/2 



P+ 



[dx 



3e 



-ip+ £ x i 



'*s(aO, (36) 



where we suppressed, for brevity, the light-like vector in 
the arguments of the fields, i.e. u+(zx) = t/lijizi), etc. 



Making use of ( 13 ) and the explicit expression for the 
proton state in ( 22 1 one finds after a short calculation 



(x 1 ,x 2 ,x 3 ) = ® 3 (x 1 ,x 2 ,x 3 ;li ) , 



(37) 



i.e. the function (j>{xi) which enters the definition (23) of 
the three-quark LCWF is nothing but the leading-twist 
nucleon DA. 

The DA $3(2;; /1) can be expanded in eigenfunctions of 
the one- loop evolution kernel Pk{x) such that the coeffi- 



cients Ck(fJ-) have autonomous scale dependence: 

00 / a (a) \ lk ^° 
$ 3 (x;ti) = 120x lX2 x 3 ( c^y) c k (n )P k (x). 

The eigenfunctions P k (x) form a specific set of homoge- 
neous polynomials of three variables which are orthogo- 
nal with respect to the conformal scalar product [3"4"] : 

120 j [dx] 3 xix 2 x 3 P k (x)Pj(x) = v k 5 kj , (39) 

where the coefficients v k depend on the normaliza- 
tion convention for the eigenfunctions P k {x). One can 
show that all eigenfunctions have definite parity under 
the interchange of the first and the third argument: 
P k (x 3 ,x 2 ,xi) = ±Pk(xi,X2,x 3 ). The first few terms in 
this expansion are 24 



$ 3 (x 1 ,x 2 ,x 3 ) = 120iAr:z:i2;2iE3 1 + a-{x x - x 3 ) 



+ b-{xi + x 3 — 2x 2 ) + 



(40) 



where we have changed the notation to Jm = cq, a — 
c\/c{), b — c 2 /cq. The corresponding anomalous dimen- 
sions are 70 = 2/3, 7 a = 20/9 and 7^ = 8/3. 

The normalization constant is determined by the 
matrix element of the corresponding local three-quark 
operator. It was calculated several times in the past using 
QCD sum rules [SBffiD]: 

f N = J [dx] 3 <S>3(x) = (5.0 ± 0.5) x 10" 3 GeV 2 . (41) 

The latest estimates for the "shape" parameters a, b from 
lattice calculations [4T1 FI2"] are in the range 



-a = 0.85 - 0.95. 



1 



b = 0.23 - 0.33. 



(42) 



These values are consistent with the light-cone sum 
rules for nucleon electromagnetic form factors |39j and 
somewhat smaller than the earlier QCD sum rule esti- 
mates [3BH3S]. 

The model used in Ref. [2TJ [22] corresponds to a = 
b = 1 at the scale /Xo = lGeV which does not contradict 
(42 1. The overall normalization constant was determined 



in Ref. [21] from the fit to parton distributions at large 
values of Bjorken x: /jv = 4.7 x 10~ 3 GeV 2 , in a remark- 
ably good agreement with Eq. (41). This agreement is 



very encouraging as a strong indication for the selfconsis- 
tence of the whole approach. Note that the coupling /jv 
is related to the normalization constant f 3 used in [2T|[22] 
as f 3 = y/2f N . 

The quark-gluon twist-4 nucleon DAs were introduced 
in [H] 
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(0\ige i ^(z 1 )u^(z 2 )[f ++ (z 4 )di(z 3 )] k \p,X) =- -m N p 2 + Nl(p) J [^] 4 e^+ 3%{x) , 

(0\ig& k *$(zi) [U+iz^u^Wdfizs) \p,X)^-\m NP 2 + Nl(p) J [dx]^** 9%{x) , 

(0^ k [f ++ (z i )ui(z 1 )Yu^(z 2 )di k (z 3 )\p, A) = - l -m N p\Ni(p) J [dx] A e^+ SfOr) , (43) 



where we changed an overall sign because of the different 
definition of the charge conjugation matrix, cf. [33] . 
The asymptotic DAs are 

$f(x,/i) = -^Ix^xsxj Af(/i) - -Af(/i) 



Sf(x,/i) 



5!xia;2X3a: 4 



Af(M) + ^Af( M ) 



>! ^1X2X30:4 Af (/i) 



(44) 



where Af are multiplicatively renormalizable couplings 



Af( A1 ) = Af( Al0 )L 7 ^, 
Af( M ) = Af( M o)i 79/(9ft) 



(45) 



with L = a s (n)/a s (fMo) and /?o = H — | n /- I n notation 
of Ref. [24] Af = ff , Af = r& fi and Af = 

Numerical values of these parameters can be estimated 
using QCD sum rules, see App. A. We obtain at the scale 
~ 1 GeV: 

Af =(2.6 ± 1.2) ■ l(T 3 GeV 2 , 
Af =(2.3 ± 0.7) ■ 1CT 3 GeV 2 , 
(0.54 ±0.2) • 10~ 3 GeV 2 , 



^3 



(46) 



where the sign convention is that the three-quark cou- 
pling /jy is positive. 

Evaluating the matrix elements in the definitions of 
DAs (43 1 using ( 13 1, (151 and explicit expressions for the 



uudg Fock states in terms of the corresponding LCWFs 
one obtains the required relations: 



- mN 
96 



2Sf (ari, x 2 , x 3 , i 4 ) + El(x 2 , xi,x 3 , x 4 ) 



m N 
' 48 



^l(x 2 ,Xx,X 3 ,X4,) + ~$l(xi,X2,X 3 ,X4) 



ilfjf ) (x 1 ,X 3 ,X 2 ,X4) 



1 



® 9 4 (xi,x 2 ,x 3 ,x 4 ) + -^l(x 2 ,x 1 ,x 3 ,x 4 ) .(47) 



48 

Note that the DA S| satisfies the symmetry relation [24] 

Zl(x 1 ,x 2 ,x 3 ,x 4 ) +E 9 4 (x 1 ,x 3 ,x 2 ,x i ) = 

= El(x 2 ,X 3 ,Xl,X4) +EI(X3,X2,X1,X4,) (48) 



which is consistent with Eq. ( 28 ) 



C. Fock state probabilities 

Our conventions correspond to the usual relativistic 
normalization of the proton state 

(p, +\P', +) = W 2P+S(P+ - P+) S 2 (p± - p'J . (49) 

The partial contribution of each Fock state is defined 
similarly, e.g. 



(p,+\p'>+)uud = 

= (2ir) 3 2p+S(p^ 



P ' + )5 2 (p ± -p> ± )P u 



(50) 



where P uu d is the probability of the three-quark state 
with zero orbital angular momentum. 

Using the definition in Eq. ( 22 ) and the ansatz in 



Eq. ( 23 ) we get after the integration over transverse mo- 



menta 



P, 



[dx] 3 

XiX 2 X 3 



\cj>(x 1 ,x 2 ,x 3 ) ] { i 



+ 2 \<t>( x 3> x 2 ,xi) + 4>{x 1 ,x 2 ,x 3 )Y 



(51) 



where p 3 = pn=3 is defined in Eq. (25 1. 

For the model specified in Eq. (38) one obtains 



15 f2 



1 



56 



(52) 



For a given value of the wave function at the origin, 
/at, the probability of the thrcc-quark valence state is 
proportional to the fourth power of the 03 parameter, 
Puud ~ flf- We fix 03 from the requirement to have 
the same probability of three-quark state as in [211 E2] . 
Namely, for f N = 5 x 10~ 3 GeV 2 and a = b = 1 one gets 



Puud 



435 
112 



JnP3 



0.17 



for 



a 3 = 0.73 GeV^ 1 . 



(53) 



(54) 



The dependence on the shape of the DA (for a,b ~ 1) 
is very weak. This property is due to an attractive fea- 



ture of the Bolz-Kroll ansatz ( 23 ) : Different terms in the 
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expansion of the DA in multiplicatively renormalizable where 0,3 = 0.73 GeV 1 , Eq. (54 1. The choice a^. = a J 



operators ( 38 1 contribute to the norm additively; there 



is no interference. For the general case one obtains 
5 

k 



(55) 



where are the expansion coefficients corresponding 
to the eigenfunctions P k with positive (negative) par- 
ity with respect to the permutation x 4 -f-> 2:3: Each state 
with positive parity contributes with extra factor three. 

The probabilities of the four-parton states with an ex- 
tra gluon with negative (positive) helicity are given by 



<-dg i 



P, 



uudg^ 



2 P4 



Ap 4 



[dx] 4 
w(x) 
[dx] 4 
w(x) 



<fi g (x)(2-Vi 2 )<j> g (x) 
{(^P+V^fKx) 



(56) 

where we use a shorthand notation 



w(x) = X1X2X3X4 



(57) 



and 



ip g (x) = i/j g 1) (x 1 ,x 2 ,x 3 ,x 4 ) - ip g 2 \x 3 ,Xi,x 2 ,x 4 ) . (58) 

Here and below V12, 7^23 , etc. are quark permutation 
operators, e.g. T'i2ip g {xi,X2 ) x 3 ,x A ) = ^ g (x 2 ,x 1 ,x 3 ,x 4 ). 



We also assumed that the functions 
real. 



and ip 



(1,2) 



arc 



With the help of Eqs. (47 1 one can rewrite these ex- 
pressions in terms of the nucleon DAs, H|, and $|. 



Using asymptotic DAs specified in (44) and the value 
a s = 0.5 (at the scale 1 GeV) one obtains for the central 
values of the couplings in Eq. ( 46 1 



a 3 corresponds to the same spread in transverse plane as 
for the three-quark wave function. These numbers are of 
the right order of magnitude, which is encouraging. 

For the general case, the DAs S|, VE^ and $4 can be 
expanded in contributions of multiplicatively renormal- 
izable operators as follows [24] 



El(x 1 ,x 2 ,X3,x i ;fi) = (t) Q (x)^2cf(fj,) P^(x) , 



^l{x 1 ,x 2 ,x 3 ,x 4 ;fj,) ± $1(x 3 ,Xi,X2,x 4 ;(ji) 



<Mx)Vc±( M )P±(x),(60) 



where <t>o{x) = \%\x\x 2 x 3 x\ and P k ~ J ' ± (x 4 , x 2 , x 3 , x 4 ) are 
orthogonal polynomials which we assume normalized as 



u~8 kr = / [dx] 4 4, (x) Pjf(x) (2 + V l2 ) P?(x) , 
v±8 kr = f [dx} 4 4> {x) If (x) (2 ± V 23 ) Pr [x) ■ (61) 



Inserting (47) and (60) into (56) one finds after some 
algebra 



_ 105?n^p 4 \ ^ S| 3| 

^uudg± — ^2 2-^i V k\ c k\ 



P 



uudg^ 



8g 2 

105m^rP4 



(62) 



Similar to the three-quark case, each multiplicatively 
renormalizable contribution to the DA generates an ad- 
ditive contribution to the state probability; there is no 
interference. 



uudg^ 



35 



105 



P4(A? 



N 



: m N p 4 



g\2 




(A|) 2 + (Ag) 2 



0.37 




(59) 



IV. PARTON DENSITIES 

The definitions of quark and gluon parton densities can 
be found e.g. in the review [48] . Translating them into 
the two-component spinor notation we obtain for quark 
and gluon distributions 



q{x) = 2 / ^ ei " (pn) ^^(-^ n )«i(^ n )+«t(-^ n )3+(^ n )if)' 

Aq{x) = \ I ^ eiXZ(Pn) (P'+l5+(-^^)9+(^")-9+(-^^)d(^")|P>+). 
Sq(x) = \ J ^e*'W<p,Mql(-lzn)&(\zn)\p,-), (63) 
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X9{X) = 2pn/ £^^(p|^ + (-izn)/« + (izn) + /« + (-^n)/« + (izn)b), 

^A 5 (x) = 2^/ ge^")( P ,+|/; + (-^n)/; + (^n)-/; + (-^n)/^ + (^n)b,+), (64) 

I 



respectively. Here q(x), g(x) are unpolarized and 
Aq(x), Ag(x) polarized densities, and Sq(x) is the quark 
transversity. For the unpolarized distributions the aver- 
age over the proton polarizations is assumed. 

The quark parton distributions for each flavor q = u,d 
receive contributions from the three-quark 3q = uud Fock 
state and also from the 3qg = uudg states with both 



gluon helicities: 

= Q3q(x) + q 3ggi (x) + q 3qg ^(x) , ( 65 ) 

and similar for Aq(x) and 5q(x). The three-quark con- 
tributions are: 



( U 3q(x) 

\Au 3q (x) 

( d 3q (x) 
\Ad 3q (x) 

and 



Pzfli f 1 dx 2 dx 3 



96x J Q x 2 x 3 

N 



96x 



5(1 -x-x 2 - x 3 )^<j) 2 (x, x 2 ,x 3 ) ± 4> 2 (x 2 , x, x 3 ) + <j)(x, x 3 , x 2 ) + cp(x 2 , x 3 , x)J j, 

(66) 



5u 3q (x) 



j dxidx^ _ x _ Xi _ X2 )f ( j ) 2 (x 1 ,x 2 ,x) ± \\<j>(xux,x 2 ) + <t>{x 2 ,x,x 1 )] }, 
Jo xix 2 L 2 L J J 

5(1 - x - x 2 - x 3 ) <j>(x, x 2 ,x 3 ) + <j)(x 3 ,x 2 ,x) <t>(x, x 3 , x 2 ) + <j>(x 2 , x 3 , x) 



Psf% f 1 dx 2 dx 3 



96x J x 2 x 3 



Sd 3q (x) = -^P^ [ dxidx2 s(l -x-x-i- x 2 )(j>(x 1 ,x 2 ,x)(f>(x 2 ,xi,x) . 
96x J x x x 2 



For the three-quark-gluon contributions we obtain: 



( U 3q 9l (x) 

\Au 3qgx (x) 

( d 3q gi (x) 

\Ad 3qgi (x) 

( u 3qgt (x) 
\Au 3q3f (x) 



( d 3qgt (x) 
\Ad 3qgf (x) 



= 2p 4 



[dx}4 
w(x) 



[5(x - an) + 5(x - x 2 )]4> g (x) (2 - V 12 ) </> g (x) , 



2p4 J W) 5[x ~ X3)Mx) (2 " Vl2) Mx) 

ix]i 
Ax) 



2 P 4 I Sy{2[<5(^ -x 2 )±5(x- Xl )] [(^(x) + V 23 ^f\x)Y - ^(x)V 23 ^ 2 \x) 



± 



[5(x - x x ) + 5(x - x 3 )]i) g (x) (2 - Via) M x )}> 



2 P4 



[ ^{s(x x 2 ) i, g (x) (2 - V l3 ) 4> g (x) ± 25(x x 3 ) [(V«(x) + V 23 ^ 2 \x)) ' 
-^(x)V 23 ^f\x)]}, 



and 



6d 3qg ,(x) = 6d 3qgl (x) = -2 P4 J [ ^5(x-x 2 )(v 23 <t> g (x))(2-V 13 ) 



M x ), 



5u 3qg r(x)=5u 3qg ,(x) = 2p 4 I ^J(x - x 2 ) [^\x){2 - V 12 )<t> g (x) + (v 23 ^f\x)) (l + Pi 2 )<t> g (xj\ . (69) 



(67) 



(68) 



•)(x) 

Finally, for the gluon parton distributions we get 



( 9(x) 
\Ag(x) 



4P4 / [dx] ^ X x) XA) {{^\x) + V 23 ^(x)) 2 - ^ 2 \x)V 23 ^\x) + \^ g (x){2 V 13 )^ g (x) 



± 



4> g (x)(2-V 12 )4> g {xj\}. 



(70) 



r 



For simple models of the wave functions the integrations plicitly. In particular using the three-quark wave func- 
over parton momentum fractions can be carried out ex- 
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which corresponds to the choice 



tion from Ref. 
a = b = 1 in the nucleon DA ( |40[ ), one obtains 

1960 



U3q(x) 
d3q(x) 



= P; 



, tyou , 3 r b., . u M , 2 i 
(i-^li + sa-^ + ya-x) 2 }, 



140 



12 
12. 

. , . _ 5600 „ . 3 r 21 .9 M . 2 i 
A<x 3g (z) = P 3? — x(l~x) 3 [l- — (l-z) + — (1-x) 2 }, 

140 r 9 ~i 

Ad 3g (z) = -P 3g ~ x{l-x) 3 {l+3(l-x) + ~(l-x) 2 j, 



fe 3g 0r) - P 39 ^ r ^(l-a ; ) 3 {l-5(l-a:) + ^(l-a ; ) 2 }, 



^^^(l-xfji+sci-^+la-x) 2 }. 

(71) 
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These expressions coincide with the corresponding ones 
in Ref. [2TJ [22]. For the three-quark-gluon contribu- 
tions, taking into account Eqs. (47) and using asymptotic 
D As ( 44 ) we arrive at 



d 3q gt(x) = -u 3qgt (x) = 56P 3qgt x(l-x) 6 , (72) 



Ad 3qg ^(x) 
Ad 3qg ^ {x) 
Au 3qg t(x) 



-Au 3qgi (x) = d 3q9i {x) . 



1 - ) d 3qg t(x), 



P u 3qg^ (x) 



(73) 



Su 3qg t(x) = fc 39S 4(a;) 

= 84 ^ - l/3v^) y/P 3qg ,P 3m x(l - xf 
5d 3qgf (x) = <W 3(?ff i(a;) 



and 



i(;(a:) = 168(P 3gg t+P 3m )a; 3 (l-x) 5 , 



(74) 



zA^x) = 168(P 39ffT -P 3m )a ; 3 (l- a; ) 5 . (75) 
where we used the notation 

= n?i¥n^2 = °- 052 ± °- 030 ■ ( 76 ) 
(^2) + 

It is easy to check that the Soffer inequality [SH [53] 

g(s) + Aq{x) > 2\Sq(x)\ (77) 
is fulfilled for arbitrary values of the parameters. 



Note that our result for the the large x behavior of 
quark parton distributions due to contribution of the 
quark-gluon Fock states differs from that in (3TJ: (1— a;) 6 
vs. (1 — x) 7 . 

For the numerical analysis we accept the same three- 
quark wave function as in Refs. [211 122] . corresponding 
to the probability of the valence state P 3q = 0.17 ( p33| , 
and fix the remaining parameters of the the quark-gluon 
wave functions from the requirement that the resulting 
parton distributions are in reasonable agreement with the 
existing parameterizations at large x, see Fig. [T] 

The unpolarized distributions are only sensitive to the 
total probability to find an extra gluon. We choose 



P 3 qg — P 3 qgf + P^qg^ — 0-33 ■ 



(78) 



For the central values of the QCD sum rule estimates 



for the wave functions at the origin, Eq. ( 46 ) , this value 



can be obtained assuming that the quark-gluon state is 
slightly more compact in transverse space as compared 
to the valence three-quark configuration: 



0.9 a 3 , 



(79) 



which is reasonable. 

The ratio Ag(x)/g{x) is determined in our simple 
model by the ratio of the probabilities to find a gluon 
with helicity aligned an anti-aligned with that of the pro- 
ton. In the rest of this work we take 



P.. 



(A? 



9^2 



P 



uudg^ 



3 (Af)* 



(Af) 2 



0.6 (0.8 ±0.2) (80) 



where the number in parenthesis is the QCD sum rule 
prediction, Eq. (46). The polarized quark distributions 



Au(x) and Ad (x) also involve another ratio of the cou- 
plings, cf. Eq. (76), which is, however, small according to 



our estimates. The corresponding contributions to Au(x) 
and Ad(x) are below 5%. 

The results for the transversity distributions Su(x), 
Sd(x) are shown in Fig. [2j These distributions are only 
very weakly constrained by the experiment, see e.g. the 
discussion in Refs. | 5~iTi r 56] . Our results are generally sim- 
ilar to the other existing model predictions, see Ref. [57] 
for a review and the corresponding references. 

We remind that in this work we try to keep the model 
as simple as possible, restricting ourselves to contribu- 
tions of the states with total zero angular momentum 
and the simplest, asymptotic shape of the four-particle 
quark-gluon proton distribution amplitude. It is seen 
that this simple approximation captures main features 
of parton distributions at large x surprisingly well, al- 
though more sophisticated models are certainly needed 
for a quantitative description. 
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FIG. 1: Quark and gluon parton distributions. The black curves correspond to the existing parametrizations: GRV 26 (short 
dashes), DSSV gS] (long dashes) and LSSTO [50] (dash-dotted) at the scale fi 2 = lGeV 2 [51!. The solid blue curve is our 
model prediction taking into account the contributions of the valence three-quark state and the state involving one additional 
gluon. The contribution of the valence state alone is shown by dots for comparison. 



V. TWIST-3 OBSERVABLES 



A. Quark-antiquark-gluon correlation functions 



A description of twist-three observables in the frame- 
work of collinear factorization involves quark-antiquark- 
gluon correlation functions which are defined as matrix 
elements of nonlocal (light-ray) three-particle operators. 
In the literature there exists apparently no "standard" 
definition of such operators, and also no standard no- 
tation. One of the usual choices jUj is to consider the 



operators 

S±(zi,z 2> z 3 ) = ^q(z 1 )[iF l , + (z2) ± F t _ l+ (z 2 )j 5 ]j+q(z 3 ) 

(81) 

and define the twist-three correlations functions as 
the matrix elements 

(p,s\S* (z 1 ,z 2 ,z 3 )\p,s) = 4mNi(jm)[s M (pn) - p M (sri)] 
x J Vxe lpn ^ z ^Df(x t ), (82) 

where s M is the proton spin vector which we assume to 
be normalized as s 2 = — 1. This formulation is often 
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1 . 




. 00 



0.8 1.0 



FIG. 2: The quark transversity distribution 5q(x). The solid blue curve is our model prediction taking into account the 
contributions of the valence three-quark state and the state involving one additional gluon. The contribution of the valence 
state alone is shown by dots for comparison. The Soffer bound ( 77 1 is indicated by the (magenta) dashed curve. 



used e.g. in the studies of the nucleon structure function 
52 (z, Q 2 )- 

Here and below the integration measure T>x is defined 

as 



Vx = 



dx\dx 2 dx 3 S(y Xi) 



(83) 



The difference to [da;] 3 (20) is that the momentum frac- 
tions sum up to zero. 

A subtlety in using this definition is that the twist- 
three and twist-four contributions in are not sepa- 
rated on the operator level. It can be more convenient to 
forgo the explicit Lorentz covariance and restrict oneself 
to transverse spin polarizations (st • n) = introducing 
another set of operators [17] : 

6 ± (z) = 2is T [§+(z 1 ,z 2 ,z 3 )±§-(z 3 ,z 2 ,z 1 )]. (84) 

The operators & ± (z) are even (odd) with respect to the 
charge conjugation. One can show that [47] 

{e ± (z)) i = ±e ± (z) 

so that the C-even ("plus") and C-odd ("minus") 6- 
operators are hermitian and antihcrmitian, respectively. 

The corresponding matrix elements define the C-even 
and the C-odd twist-three correlations functions 

(p,s T |© ± (^)b,ST) =2(pn) 2 J Vxe- l(pn) ^ XkZk e ± (x), 

(85) 

which are related to the D ± -functions introduced above 
as 

8m N D + (x 1 ,x 2 ,x 3 ,) = & + {x 1 ,x 2 ,x 3 ) - &~{x 1 ,x 2 ,x 3 ) , 
8m N D~ (x 1 ,x 2 ,x 3 ) = & + (x 3 ,x 2 ,xi) + &~(x a ,x 2 ,Xi) . 

(86) 

Note that we use the same notation & ± for the operators 
and the matrix elements, which hopefully will not lead 
to a confusion. 



The helicity structure of the twist-three correlation 
functions can be made explicit going over to the two- 
component spinor notation. One obtains 



@ ± (z) = -ig sQ ± (z)-sQ ± (z) 
where s — s 1 + is 2 , s = s 1 — is 2 , and 



(87) 



Q ± (z) = qi(z 1 )f ++ (z 2 )qi(z 3 )±ql(z 3 )f ++ (z 2 )ql(z 1 ), 



Q*{z) - q[(z 1 )f ++ (z 2 )q[(z 3 )±q l + (z 3 )f ++ (z 2 )qX(z 1 ) 



The nucleon state with a transverse polarization can be 
expressed in terms of the helicity states \p, ±) as 



\p,st) 



1 

T2 



s \p,-) 



Taking into account that the operators Q ± (z) increase 
and Q ± (z) decrease helicity, it follows that 

(p, stI©* (z)\p, s t ) = 

= ig[(p,+\Q ± (z)\P,-)-(P,-\Q ± (z)\P,+)]-(m 

It is easy to see that Q (z) — ±[Q ± (z)]^, so that the 
two matrix elements on the r.h.s. of Eq. (89) are related 



and one does not need to consider the operators with a 
"tilde" explicitly. 
We define the distributions Q ± (x) as 

(p,+|C ± (2)|p,-) =-2i{pn) 2 J Vxe- l(pn ^ XkZk Q ± (x). 

(90) 

The P-parity implies (cf. [47] ) 



and finally 



(Q ± (-x)Y =±Q ± (x) 
G ± (x) = -gQ ± (x). 



(91) 
(92) 




In the light-cone formalism, twist-three correlation 
functions are generated by the interference of Fock states 
with different particle content as illustrated in Fig. [3] 
The contributions shown schematically in Fig. |3^,,b cor- 
respond to the interference of the three-quark and three- 
quark-gluon wave functions, whereas the one in Fig. [3p 
stands for the interference of the three-quark-gluon state 
with the one containing an extra quark-antiquark pair. 
The latter term contributes to a different kinematic re- 
gion in momentum fractions compared to the first two 
terms and is missing to our accuracy. 

Explicit expressions for the three-quark and three- 
quark-gluon Fock states for the nucleon with positive 



helicity are given in Sec. Ill The corresponding states 
for the nucleon of negative helicity are given by the 
same expressions, Eqs. (22 1, (271, (31), where helic- 
ities of creation operators have to be flipped. The 
wave functions of the three-quark-gluon states of the nu- 
cleon with positive and negative helicity are the same, 
1*1234] (-) = 1*1234] (+)' whereas for the valence three- 
quark state there is an overall sign difference: [^123] (-) = 



,(0) 

12; 
3se( 
defined as 



1*123] (+)• All matrix elements in question can be ex- 

t( 



pressed in terms of two correlation functions Q q (x) 



uud(P, +\<l+( z 3).f++( z 2)q+(zi)\p, -)uudgt 

= -2ip\ J Vxe- tp +^ x ^ Q\{x), 

uud(P, +\Q+{ Z l)f++( z 2)q+{Z3)\p, ->W ff T 

= -2ip\ J Vxe- tp +^ x ^Q\{x), (93) 

where the subscript q = u,d stands for quark flavor. In 
particular 



q£(x) ± sic-*) + Ql(-x) ± Q\ 



(94) 



Using the ansatz for the LCWFs in Eqs. (23), (32 1 one 



can represent (x) as convolution integrals of the dis- 
tribution amplitudes. We obtain 



Qd( X ) =\ A 9{- X 1, X 2-> X Z)~ [ ^^^(i + ^l-Cl-6) 0(6,-^1, 6) + 0(6,-^1,6) i>g(€l,X3,&,X2), 

& X\ J qi t,2 



1 



Qt(x) =-A9(- Xl ,x 2 ,x 3 ) — 



1 fd£id& 



x\ J 6 6 



8(i + x!-ti- 6)0(6 , -xi , 6) [v^ (6 , x 3 , 6 , x 2 ) - v4 2) (6 ,b,x 3 ,x 2 ) 



Q\{x) =Ae( Xl ,x 2 ,-x 3 )- I ^^6(1 + x 3 -6 -6)0(6,6, -a*) [^^(6,6^1^2) + ^ 2) (6,^i, 6,^2) 
fit (*) = - A6{ Xl ,x 2 , -x 3 )± J ^ ^5(1 + x 3 - 6 - 6){0(-*3, 6,6) ^(n, 6, 6,^2) + \l>f> (n, 6, 6, x 2 ) 
(-3:3,6,6) + 0(6,6, -^3) ip g (xi, 6 ,6, x 2 ) - 2^(6,6,2:1,^2) }, (95) 



where it is implicitly assumed that £1 + #2 + £3 = 0, the 
Heaviside step- function with several arguments is defined 
as 9(a, b, c) = 0{a)8{b)6(c), and 



(96) 



The distribution ip g is defined in Eq. (58 1 . The QCD sum 

ifrip and 



rule result Af <C Af 2 



(46) implies that tp 



(i) 



as a consequence both helicity down functions d {x) 
are suppressed in comparison with the helicity up func- 
tions Q} u( i( x )- Note that we use a symmetric notation 
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FIG. 4: Support properties of twist-three correlation func- 
tions in barycentric coordinates. For the explanation of dif- 
ferent regions see text. 



with quark, antiquark and gluon momentum fractions 
are treated equally so that the momentum conservation 
condition is x\ + x 2 + x 3 = 0. Support properties of the 
correlation functions (58) can most easily be shown going 
over to barycentric coordinates |47j as shown in Fig. [4] 

x = x\e\ + x 2 e 2 + x 3 e 3 = X\E\ + x 2 E 2 . 

Three-parton correlation functions, in general, "live" in- 
side a hexagon-shaped area which can further be decom- 
posed in six different regions (triangles). The triangles 
labeled (12) + 3~, 2 + (13)~, etc., correspond to different 
subprocesses at the parton level [55]; For each parton 
k = 1,2,3 "plus" stays for emission (xu > 0) and "mi- 
nus" for absorption (xj. < 0). Alternatively, one may 
think of "plus" and "minus" labels as indicating whether 
the corresponding parton appears in the direct or the fi- 
nal amplitude in the cut diagram, cf. Fig. [3| It is impor- 
tant that different regions do not have autonomous scale 
dependence; they "talk" to each other and get mixed un- 
der the evolution, see Ref. [17J for a detailed discussion. 

Our model predictions for the correlation f unc tions 
Qd(x) and — Qu(x) (note opposite sign), Eq 
shown in Fig. [3] and Fig. |6j respectively. Bot 



(|90J), are 
distri- 
butions are symmetric with respect to the center of the 
hexagon: Qg(xi,x 2 ,x 3 ) = Qq(-xi,-x 2 ,-x 3 ), which is 
a consequence of P-parity, cf. Eq. (91 ). Each of the four 



terms Qj (±2;) in Eq. (94) is confined to a different 
"triangle" and, hence, has a different partonic interpre- 
tation: 



Qj(z) 



(12)+3- 

(23)+r 



-x) 
-x) 



3+(12)- 
l+(23)- 



(97) 



The larger contributions, e.g. in the (12) + 3~ region, 
correspond to (valence) quark emission with momentum 
fraction X\ > and subsequent absorption with momen- 
tum fraction — x 3 > X\ > 0, accompanied with gluon 
emission with momentum fraction x 2 > 0. The smaller 
contributions, e.g. in the l + (23)~ region, differ from the 
above in that the gluon with momentum fraction — x 2 > 
is absorbed and thus x\ > —x 3 > 0. Note that there is no 



symmetry between gluon emission and absorption, which 
may be somewhat counterintuitive. 

The dominant, gluon emission contribution to the uGu 
correlation function Q\{x) is roughly factor two larger 
compared to the dGd distribution, Q/Ax), and has the 
opposite sign. The contributions of gluon absorption, 
Qu(x) and Q*,(x), have the same sign for u- and d-quarks, 
and are much smaller compared to gluon emission. 

Our model correlation functions vanish in the 2 + (13)~ 
and (13) + 2~ regions. This property is an artefact of ne- 
glecting contributions of the type shown in Fig. [3J; which 
are formally higher order in the Fock expansion. These 
contributions can be estimated using a model for the five- 
parton qqq{qq) state from Ref. [21] and turn out to be 
considerably smaller than the ones considered here. 

The "minus" correlation functions Q„(x) and Qj(x) 
are obtained from the "plus" ones by changing the sign 
of the contributions in the (12) + 3~ and l + (23)~ regions, 
so we do not show them separately. 



B. The structure function g2(x,Q 2 ) 

The structure function g 2 (xB, Q 2 ) is given by the sum 
of the Wandzura-Wilczek (WW) and genuine twist-3 con- 
tributions 

92(x B , Q 2 ) - g? w (x B , Q 2 ) + gl w - 3 (x B , Q 2 ) . (98) 
The WW contribution reads 



g? w (x B ,Q') = - gi (x Bl Q 2 ) + 



where 



1 dy 



5l (y,Q 2 ), (99) 



gi{x B ,Q 2 ) = \ 4 [&q{x B ,Q 2 ) + Aq(-x B , Q 2 )] . 

(100) 

The twist-3 contribution g%" 3 (x B ,Q 2 ) can be written 

as 



gT-*(x B ,Q') = 
1 2 f 1 d£ 



Aq T ^,Q 2 ) + Aq T (-C,Q 2 )[W) 



where A(?t(£) is defined in terms of the D + - function 
introduced in Eq. ( 82 ) : 



A?t(£)=4 j VxD+(x 



t j_ \ ^ + x 3 )-s^- Xl y 

dx 3 [ x\ + x 3 

(102) 

As above, the subscript q refers to the contribution of a 
given quark flavor. In terms of the Q q (x)-functions 
one obtains 



D+(x) - 



4tojv 



[Qj(s) + Qi(- 



(103) 
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FIG. 5: The quark-antiquark-gluon twist-three correlation function Q^(x) at the reference scale 1 GeV. 
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FIG. 6: The quark-antiquark-gluon twist-three correlation function —Qt(x) at the reference scale 1 GeV. 



To avoid confusion, in this section we use the notation 
xb for the Bjorken variable, whereas x is reserved for the 
set of parton momentum fractions x = {x\, £2, 2:3 }■ 

Under a plausible assumption that the spin-dependent 
part of the forward Compton amplitude satisfies a dis- 
persion relation without subtractions, the integral of 
g2(xB,Q 2 ) and, hence, of g%"~ (xb,Q 2 ) vanishes [S"§] 



f dx B g t 2 w - 3 (xB,Q 2 ) = 0. 
Jo 



(104) 



This statement is known as the Burkhardt-Cottingham 
(BC) sum rule. 



Using Eqs. (95) one finds that in our model D+(x) 



is nonzero only when x± > and X3 < 0. This, 
in turn, implies that A^t(C) vanishes for £ < (i.e. 
there is no antiquark contribution). As a consequence, 
in our model g%"~ 3 {xb,Q 2 ) satisfies in addition to 
Burkhardt-Cottingham (104) also the Efremov-Leader- 
Teryaev (ELT) sum rule [§Up 



dx B x B g t 2 w - 3 (xB,Q 2 ) = 0. 



(105) 
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FIG. 7: Upper panels: Experimental results on the proton (left) and neutron (right) structure function qi{xb^ Q 2 ) compared to 
our model calculation at the scale Q 2 = 1 GeV 2 . Lower panels: The twist-three contributions xgl wlst ~ 3 (xB,Q 2 ) for the proton 
(left) and neutron (right) compared to the analysis in Ref. |9J (shaded areas). Our model predictions at the scale Q 2 — 1 GeV 2 
and Q 2 — 10 GeV 2 are shown by the black solid and dashed red curves, respectively. The predictions at 10 GeV 2 obtained using an 
approximate evolution equation from Refs. [6 68J are shown by the red dotted curves for comparison. 



For the second moments one obtains 

d 2 , p = 3 / dx B x 2 B g t 2 w ~ :i (x B ) 
Jo 

5 a, l\a A 5a + b\ .( a + 56\i 

= 3^M A H 1+ ^J +Af ( 1+ ^H]' 

d,2,n = 3 / dx B x B gl w ~ 3 (x B ) 
Jo 

5 ap l\aA. b — 5a\ _/ a — 5by 

(106) 

The corresponding numerical values are, at the scale 
1 GeV: 

d 2 , P = 0.0016, d 2 , n = -0.00072. (107) 

Both numbers compare very well to the lattice QCD [65] , 
QCD sum rules [5l[66] and chiral quark soliton model [67] 
calculations. The negative value of d% for the neutron 
(in all models) is in conflict, however, with the existing 



experimental average: 

d^p = °- 0032 ± °- 0017 EH , 

d^P = 0.0062 ± 0.0028 [§3] . (108) 

Further, a straightforward calculation gives 

9^ P ~ 3 M = 0-0436772 ( In x B + x B + 1/2b|) 

+ x% (l.57357 - 5.94918x B + 6.74412x| - 2.19114z|) , 

S^n 8 M = 0.0655158(111^ + x B + l/2x 2 B ^j 

+ x% (0.130996 - 1.121015s + 2.31342x| - 1.20598x|) 

(109) 

(at the reference scale 1 GeV) for the proton and neutron, 
respectively. 

Our results for the full structure function g2{x B ,Q 2 ) 
are compared to the experimental data [62H64] in Fig. [7] 
(upper panels) and, separately, for the twist-three con- 
tribution g%"~ 3 (xB) to the analysis in Ref. [S] (lower 
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panels). The twist-three contributions are shown at the 
model scale Q 2 = 1 GeV 2 and after the evolution to a 
higher scale Q 2 = 10 GeV 2 . The scale dependence was 
calculated in two ways: using exact (one-loop) evolution 
equations for the relevant quark-antiquark-gluon correla- 
tion functions from Ref. [47] (dashed curves), and using 
the much simpler evolution equation from Refs. [51 [55] 
which is based on the large- N c and large-x^ approxima- 
tion and only involves the g%"~ 3 (xB) structure function 
itself (dotted curves). Since we are interested primarily 
in the large xb region, we used flavor-nonsinglet evolu- 
tion equations which are simpler. The results of both 
approaches almost coincide within the line thickness. A 
good accuracy of this approximation was expected but 
has never been checked in a dynamical model calculation. 
Note that effects of the evolution are generally significant 
because of large anomalous dimensions of twist-three op- 
erators, and have to be taken into account in the analysis 
of the experimental data. 

As seen from Fig. [7J the twist-three contribution to 
the structure function g 2 {xB,Q 2 ) at large Xb proves to 
be positive for the proton and negative for the neutron. 
This prediction can be traced to the relative signs of 
the three-quark-gluon couplings and is largely model- 
independent. It is in agreement with Ref. [9]. In the 
intermediate region 0.2 < xb < 0.5 the twist-three con- 
tribution gl W p~ 3 (xB) changes sign and becomes negative 
in our calculation, whereas it remains positive according 
to the data analysis in Ref. [9] . This difference may well 
be due to contributions of higher Fock states, with two or 
more gluons, and probably also partially remedied by us- 
ing a more sophisticated model for the three-quark-gluon 
wave function. A detailed analysis would be interesting 
but goes beyond the tasks of this paper. Another is- 
sue is that in [5] the ELT sum rule is strongly violated, 
which suggests the existence of a large positive flavor- 
singlet contribution at xb ~ 0.1 due to gluons or sea 
quark-antiquark pairs. Such contributions are related to 
the twist-three thrcc-gluon correlation functions and are 
missing in our present framework. 



are expressed in terms of Q q -functions as follows: 



C. Single spin asymmetries 



The quark-antiquark-gluon correlation functions con- 
sidered in this work are precisely those responsible for 
transverse single spin asymmetries (SSA) observed in 
different hadronic reactions, if described in the frame- 
work of collinear factorization 10 20 . The distributions 
TqF q (x), ATqFq(x) introduced in this context in Ref. |47j 



TqFq(xi,X 2l X 3 ) 
1 

4 



(i+v 13 )e+(x) + {i-v 13 )&-(x) 

Ql(x 3 ,x 2 ,xi) + Ql(-x 1 ,-X2,-x 3 ) 



_9 
2 

+ Q,\{xx, x 2 , x 3 ) + Q^(-x 3 , -x 2 , -xi) 



ATq Fq (xi,X 2 ,X 3 ) 

1 

~4 



(i-Vi 3 )e + (x) + (i+v 13 )e- (x) 

Q t q(x 3 ,X 2 ,X 1 ) - Ql(-X 1 ,-X 2) -X a ) 

- Q q (xi,x 2 , x 3 ) + Q q (-x 3 , -x 2 , -xi) • (110) 
A common notation [69J is to show quark momenta only: 

T q , F {x,x') = Tq Fq (-x',x' - X,X) , 

TA q ,F(x,x') = AT qFq (-x',x' - X,x) , (111) 

Written in this way, the distributions are symmetric (an- 
tisymmetric) functions of the arguments: T q .F{x,x') = 
T q>F (x',x) and AT q , F (x,x') — —AT qi F{x',x). A yet an- 
other notation for the same functions in a different nor- 
malization is used in the recent analysis in Ref. |71j : 



G q F (x,x') 



m N 
2 



TqFq ( X , X X, x) , 



G g F (x,x') = ATq Fq (-x',x' -x,x). (112) 

m N 

In the framework of collinear factorization, SSAs orig- 
inate from imaginary (pole) parts of propagators in the 
hard coefficient functions. In the leading order, taking 
a pole part enforces vanishing of one of the momentum 
fractions in the twist-3 parton distribution, and are clas- 
sified as soft gluon pole (SGP) or soft fermion pole (SFP), 
depending on which momentum is put to zero, respec- 
tively. Such "pole" contributions are therefore considered 
to be main source of the observed asymmetries and can be 
estimated from the available experimental data [701 I71j . 

Since our approximation for the nucleon wave function 
does not contain antiquarks, the T qFq , ATqFq distribu- 
tions are nonzero in the (23) + l~ and (12)~3 + regions 
only, cf. Fig. [4] Moreover, both distributions vanish at 
the boundaries of parton regions where one of the mo- 
mentum fractions goes to zero, and, hence, both SGP 
and SFP terms vanish as well. This property is an ob- 
vious artefact of the truncation of the Fock expansion 
to a few lowest components: The LCWF of each Fock 
state vanishes whenever momentum fraction of any par- 
ton goes to zero and the same property holds true for the 
correlation functions. Our model for the gluon distribu- 
tion xg(x) in Fig. [I] vanishes at x — > for the very same 
reason. 

For the leading-twist parton distributions, a possible 
way out is to assume the valence-type input at a certain 
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FIG. 8: The quark-antiquark-gluon twist-three correlation function —T^ Fd (x) at the reference scale /i 2 — 1 GeV 2 (left) and 
At 2 = 10 GeV 2 (right). 



low scale, and construct realistic dynamical models by 
applying QCD evolution equations that include multiple 
soft gluon radiation. This approach was suggested by 
Gluck, Reya and Vogt (GRV) HFJ [26] and proved to be 
very successful phenomcnologically. Exploiting the same 
idea for the twist-three distributions suggests itself. 
It is easy to see that both the SGP and SFP contribu- 



tions reappear once QCD evolution is taken into account. 
The full one-loop evolution equation for the functions 
TqF q , ATqF q is rather cumbersome and can be found 
in [47] . For our present purposes the flavor-nonsinglet 
evolution equation is sufficient. Restricting ourselves to 
the SGP kinematics X2 — > one obtains, to the one-loop 
accuracy 



T q ,F(x,x;n) = T q , F {x,x;fj, ) + — In 



Nc 
2 



2tt ~ pi I J x £ 
1 



Pqq{z)Tq, F (£,,0 + ~~2~ 9 ' F ~ ~ "> 7^,f(^' 



2 1-Z 



1 1 "I -1 M5 

TAq,F[X,0 + ~ 2z)T q , F {x,X - - ^rT Aq ,F(x,X - - N c T q , F (x,X) j , (113) 



where it is assumed that x > 0, P qq {z) is the usual 
Altarclli-Parisi splitting function and z = x/£. Even if 
T q ,F(x, x; /ig) = 0, a non-zero SGP contribution is gener- 
ated at a higher scale /x 2 . It is given by a certain integral 
of TqF q , ATqFq away from the line xi = 0, and involves 
large quark momentum fractions only, £ > x. (For a de- 
tailed discussion of integration regions in Eq. (113) see 
Ref. @7].) 

One difficulty in following the GRV approach is that 
the initial condition for the evolution has to be taken at 
a very low scale ^q RV — 0.25 GeV 2 [551125] whereas our 
model is formulated at /i 2 = 1 GeV 2 . The advantage 
of using the higher scale is that we have been able to 
use QCD perturbation theory and operator product ex- 
pansion to get some insight in the structure of the lowest 
Fock states, but the price to pay is that the nucleon at the 



scale 1 GeV already contains significant admixture of yet 
higher states, with several gluons and quark-antiquark 
pairs, which we do not know much about. These addi- 
tional contributions are not taken into account in this 
work, and this is the reason that we underestimate par- 
ton distributions at small x, cf. Fig. [T] 

A consistent implementation of the GRV program 
would require to give up QCD motivated models for the 
qqq and qqqg states and resort to purely phenomeno- 
logical parametrizations. We leave this study for fu- 
ture work. Instead, in what follows we show the re- 
sults corresponding to the evolution of our model twist- 
three parton distribution from 1 GeV 2 to an ad hoc scale 
/i 2 = 10 GeV 2 . This calculation should be considered as 
an illustration, since effects of the QCD evolution from 
the GRV scale Hq R v — 0.25 GeV 2 are, generally, much 
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FIG. 9: Radiatively generated SGP distributions G q F (x,x) at Q 2 — 10 GeV 2 rescaled by factor 10, shown by the solid curves, as 
compared to the phenomenological studies of spin asymmetries in high transverse momentum meson production in pp colisions [70ll71j. 



larger. 

As an example, in Fig. [8] we show the the quark- 
antiquark-gluon twist-three correlation function T^p^x) 
(with opposite sign) at the model scale /i 2 = 1 GeV 2 
(left) and after the evolution to /i 2 = 10 GeV 2 (right). 
As already mentioned above, in our model (left picture) 
this correlation function is only nonzero in the two left- 
most triangle regions corresponding to emission and sub- 
sequent absorption of the (valence) quark. The upper 
and the lower triangles corresponds to gluon emission 
and absorption, respectively. The longest diagonals of 
the hexagon, connecting diametrically opposite vertices, 
correspond to vanishing of one of the parton momen- 
tum fractions. In particular, on the horizontal diagonal 
X2 = 0, i.e. it corresponds to the SGP kinematics, and 
on the other two diagonals either x\ = or x 3 = 0, so 
they stand for the SFPs. The two triangles that come 
next to the right and include the upper (or the lower) 
edges of the hexagon, correspond to the contributions of 
the type shown in Fig. [3] where a gluon is emitted and a 
quark- antiquark pair is absorbed (or vice versa). These 
contributions are thus analogous to the so-called ERBL 
regions in off-forward parton distributions and, formally, 
are of higher order in the Fock expansion. Finally, the 
two right-most triangles correspond to the antiquark dis- 
tributions. 

Once the QCD evolution is taken into account, dif- 
ferent parton regions get mixed. In particular the gap 
between the X2 > and xi < regions gets closed and 
the SGP term appears, see Fig. [8] (right picture). The 
SFP terms are also generated, but remain very small be- 
cause the corresponding terms in the evolution equations 
are 1/Nc suppressed. 

The radiatively generated SGP distributions G q F (x,x) 
at Q 2 = 10 GeV 2 are shown in Fig. [9] and compared 
there with the results of phenomenological studies of 
spin asymmetries in high transverse momentum meson 
production in pp colisions [701 171) . Our distributions 



are of the same sign and similar shape compared to 
these studies, but about one order of magnitude smaller. 
It is plausible that much larger SGP contributions can 
be generated from the similar valence-like ansatz if the 
QCD evolution is started at a low scale of the order of 
H 2 GKV ~ 0.25 GeV 2 (2SH2B]. The SFP contributions that 
we obtain in this exercise appear to be two orders of mag- 
nitude below the estimates in Ref. [7T], albeit with the 
correct sign. It is unlikely that such large contributions 
can be obtained radiatively starting from the valence- like 
ansatz, unless one assumes the existence of antiquarks 
with large momentum fraction at low scales in the pro- 
ton WF. 



VI. CONCLUSIONS 

In this work we explored the possibility to construct 
higher-twist parton distributions in a nuclcon at some 
low reference scale from convolution integrals of the light- 
cone wave functions. 

To this end we have studied the general structure and 
introduced simple models for the four-particle nucleon 
LCWFs involving three valence quarks and a gluon with 
total orbital momentum zero, and estimated their nor- 
malization (WF at the origin) using QCD sum rules. We 
have shown that truncating the Fock expansion at this or- 
der, that is taking into account valence three-quark con- 
figuration and those with one additional gluon, provides 
one with a reasonable description of both polarized and 
unpolarized parton densities at large values of Bjorken 
variable x > 0.5. 

Using this set of LCWFs, twist-three quark-antiquark- 
gluon parton distributions have been constructed as con- 
volution integrals of qqqg and valence three-quark compo- 
nents, which enter the description of many hard reactions 
in QCD in the framework of collinear factorization. In 
particular the twist-three contribution to the polarized 
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structure function g2 (x, Q 2 ) is given by a certain integral 
of the three-particle distribution over the parton momen- 
tum fractions, and thus is a measure of its "global" prop- 
erties. Our calculation correctly reproduces the sign and 
the order of magnitude of the twist-3 term at large x, 
without free parameters. 

Transverse single spin asymmetries, on the other hand, 
are sensitive to "local" properties of the three-particle 
correlation functions in specific configurations where one 
of the momentum fractions vanishes. Since our approx- 
imation for the nucleon wave function only includes a 
few lowest Fock components, and since the LCWF of 
each Fock state vanishes whenever momentum fraction 
of any parton goes to zero, both "soft gluon pole" and 
"soft fermion pole" terms vanish at the scale where the 
model is formulated. They are, however, generated by 
QCD evolution that brings in multiple soft gluon emis- 
sion. Our results suggest that realistic dynamical models 
of the the twist-three distributions (and the pole terms) 
can be obtained following the GRV-like approach on the 
level of WFs, i.e. assuming that the nucleon state at a 
very low scale can be described in terms of a few Fock 
components, including the valence quarks, one additional 
gluon and, probably, a quark-antiquark pair, and apply- 
ing QCD evolution equations. 

An obvious problem with this strategy is that the 
starting scale has to be chosen very low, of the order 
of ^ RV -0.25 GeV 2 |23E5I, and thus the modelling 
of the wave functions necessarily becomes purely phe- 
nomenological. In spite of this, and the usual criticism of 



the application of perturbative QCD evolution equations 
at very low scales, we believe that such an approach has 
good chances to provide us with some intuition on the 
structure of higher-twist parton distributions in general, 
which is currently not available. This work is in progress 
and the results will be published elsewhere. 
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Appendix A: QCD sum rules for the quark-gluon 
wave functions at the origin 

The definitions of quark-gluon twist-4 nucleon DAs 



( 43 1 [24] can be rewritten in conventional Dirac bispinor 
notation as follows: 



ige 



ijk/ 



ige^ k (0\ u\{z x )C^u\[z 2 ) 7"Kte«(*4)b) 



-m N p 2 + ^(p) I [dx} 3 e^^ x ^mx), 



ige 



ijk. 



uUz 1 )C>y v +d>(z 2 ) M(z 3 )F*l(z 4 )\p) = -m N p 2 + ^(p) / [dx} 3 



3?(s), 



(A.l) 



where F kl v = n^F" (t a )y. The normalization of the DAs 
is determined by the matrix elements of the correspond- 
ing local operators. In what follows we estimate these 
matrix elements using the classical SVZ QCD sum rule 



approach [4"5] . 

To this end we define isospin-1/2 twist-4 quark-gluon 
operators: 



= o*5 £ 



rjk 



rjk 



(u i (x)Cftu :i (x))Yii>d l (x) - (u^x^iid^xfjYi^ix^F^x) , 
(u l {x)C l5 iiu l {x) S j 1 v iLd k {x) - (u i (x)Cj 5 id l (x) S jj^u k (x)']F^(x) , 
(u l (x)CY^(x)^d l (x) - (y(x)Cj"jid^x)^u l (x)\F%l(x). 



(A.2) 



Matrix elements of these operators sandwiched between 
vacuum and the proton state are related to the couplings 



introduced in Eq. ( 46 ) 



<0K(0)b) = -~(a| - -Ag) m N p 2 + i lb N{p) , 
(0|% 9 (0)|p) - |(A| + \i)m N p 2 + iN{p), 
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FIG. 10: Leading-order contributions to the OPE of the cor- 
relation functions in Eq. (IA.I 



The sum rules are derived for the correlation functions 
of rf k (x) with the three-quark operators [HJ US] 

m (x) = e^ fe [^(x)C77^ J ^)]757^ fc (^), 

V2 (x) = e llk [u\x)Ca^(x)] l5 a^d k (x). (A.4) 

The corresponding couplings are well known from numer- 
ous QCD sum rule calculations 
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(A.5) 



where the numbers correspond to leading-order QCD 
sum rule results at the scale 1 GeV, see e.g. [55] , 

In particular, we consider the following correlation 
functions: 



Tr 



Tr 



jslcfxe^mrifWfjiiOm) 
dW^<0|TK(z)rh(0)}|0) 



Tr 



7 5 /d 4 xe^(0|T{r,f( a; )r7 2 (0)}|0) 



pI^Kp 2 ) 

p\n(p 2 ) 



(A.6) 



The sum rules are derived from the matching of the QCD 
calculation of the invariant functions 11^ (p 2 ) at Euclidean 
p 2 r~j — 1 GeV 2 with the dispersion integral representation 
where the nucleon contribution is written explicitly: 



IL> 2 ) 



(A§-Ag/3)A! 



a "'JV 9 ■-> 

4 m z N - p 2 

1 , (Af + Af )Ai 
-m N 2 o 

D TO^ — p z 

1 2 (Af + Ag)A 2 

~c m N 2 9 
6 TO^, — 



(A.7) 



FIG. 11: The coupling (in units of GeV 2 ) (upper panel) 
and the ratio A^/^Aj) (lower panel) as a function of the 
Borel parameter M 2 for the central values of the conden- 
sates ( A.10[ |, q s (1 GeV)=0.5. The solid line corresponds to 
y^so = 1.4 GeV and the dashed line to ^/so = 1.6 GeV. 



and the contributions of higher states and the continuum 
are modelled in the usual way as the QCD spectral den- 
sity above a certain threshold, ^/so ~ 1.5 GeV, dubbed 
the interval of duality. On the QCD side, we take into ac- 
count contributions of perturbation theory and vacuum 
condensates of dimension 4 and 6 shown in Fig. |10| The 
leading-order contributions of dimension 8 vanish for all 
cases. 

Proceeding with the standard technique we derive the 
following set of sum rules: 



2(2 7 r) 4 (Af-Af/3)Aim 2 v e~^ 



45tt 12 



9tt 



2(27r) 4 (Af + A|)Aimfv e~M 
r6 



45?r 6 12 



40a 



s 2 

a , 



2(2tt) 4 (A? + Af)A 2 mfv e"^ 



a, c b 8a, Q 
—M 6 E 3 + -M 2 E 1 + — -a 2 , 
15tt d 4 3tt ' 



(A.8) 
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where M is the Borel parameter, 

n-1 



f-^ fc! VM 2 / 



fc=0 



and 



a 
b 



-{2ir) 2 (qq) = (0.55 ±0.06) GeV 3 , 
47r(a s F 2 ) = (0.47 ±0.14) GeV 4 . 



(A.9) 



(A.10) 



are the quark and gluon condensates, respectively, at the 
scale 1 GeV. 

For the numerical analysis we substitute the coupling 
Ai in the first two equations in (A. 8 by the square root 
of the "Ioffe sum rule" [4"4] 



where m 2 , = (qgaFq) / (qq) ~ 0.65 GeV 2 , and taking into 
account that Ai is negative (which is a convention). As- 
suming the "working window" in the Borel parameter 
M 2 ^1 — 2 GeV 2 and taking into account uncertainties 
in the vacuum condensates and the continuum threshold 
y/so = 1.4 — 1.6 GeV, we obtain the numbers given in 
Eq. (46) in the text. Taking into account that to a good 



-2Ai, we get §A? 



Af. This rela- 



accuracy A2 — _, u . .>> ^- . 
tion holds to the approximation considered here (leading- 
order QCD sum rules) independent on the values of vac- 
uum condensates and other parameters. It can, however, 
only be valid on a certain (low) normalization scale as 
the anomalous dimensions of the couplings are different, 
cf. ([451). 



2(27r) 4 |Ai| 2 m 2 v e"^ = M b E 3 +-M 2 E 1 



4 m, 

~3M 2 J' 
(ATI) 
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